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Abstract. In this work it is shown that there are some spatially homogeneous 
but anisotropic models (Kantowski-Sachs and Bianchi type-III), with a positive 
cosmological constant, for which the inhomogcncitics in the distribution of matter on 
the surface of the last scattering produce anisotropies (in large angular scales d > 10°) 
that do not differ from the ones produced in Friedmann-Lemaitre-Robertson- Walker 
(FLRW) models, if the density parameters are finely tuned. Namely, for adiabatic 
initial conditions, the Sachs- Wolfe effect in these anisotropic models is equal to the 

r 

one obtained for isotropic models: 4jf = ^ e + 2 J Q^dw. This result confirms the 

e 

idea developed in previous works that with the present cosmological tests we cannot 
distinguish these anisotropic models from the FLRW models, if the Hubble parameters 
along the orthogonal directions are assumed to be approximately equal. 
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1. Introduction 

The task of proving the homogeneity and isotropy of the Universe at large scales is not 
a simple one. It is generally accepted that the Universe is spatially homogeneous as a 
result of the so called Copernican principle, that is, if we assume that we live in a typical 
place, and since the isotropy of Cosmic Microwave Background Radiation (CMBR) can 
be used to prove that the Universe is locally isotropic, then (even though we cannot find a 
proper homogeneity test), we must conclude that the Universe is spatially homogeneous 
and isotropic. This conclusion reduces drastically the space of Solutions of Einstein 
equations, and the number of possible cosmological models. 

In this way we are led to the so called FLRW cosmological models to describe our 
observations. Despite the high level of isotropy, some authors have worked on spatially 
homogeneous and anisotropic models and proved that they might agree with present 
observations. For instance, if the classical tests of cosmology are applied to a simple 
Kantowski-Sachs metric and the results compared with those obtained for the Standard 
model, the observations will not be able to distinguish between these models if the 
Hubble parameters along the orthogonal directions are assumed to be approximately 
equal [1]. Following along the same lines, we made a qualitative study [2] of three 
axially Symmetrie metrics (Kantowski-Sachs, Bianchi type-I and Bianchi type-III), with 
a cosmological constant, to analyze which were physically permitted, when we assume 
them to be bound by a high degree of isotropy, that is, although our models were 
assumed anisotropic they could be considered to be almost FLRW, as far the shear 
Parameter is concerned, from the epoch of the last scattering to the present. Recall that 
one defines to be 'close' to a FLRW model when both parameters 



are almost zero, as was stated in [3]. Here, a ab is the shear tensor, H the Hubble 
Parameter, E ab and H ab the electric and magnetic part of the Weyl tensor, respectively 
(see Appendix E). From this analysis we concluded that these models are good 
candidates for the description of the observed Universe, provided that the Hubble 
parameters are approximately equal at the last scattering. In other words, the vanishing 
of the first parameter, E 2 « 0, is sufficient to assure a FLRW-like behavior. 

Historically, the detection of the CMBR has led to constrains in theoretical models 
in the field of Cosmology, and lent a hand to the Big Bang Solutions. Indeed, the 
observed level of isotropy of the CMBR, first detected by Penzias & Wilson [4] , provides 
strong evidence for the large-scale isotropy of the Universe, and is the best argument 
in favor of an isotropically expanding Universe. Later, more precise experiments 
proved that this radiation has temperature fluetuations, or anisotropies. These small 
anisotropies are thought to give rise to the observed galaxies, and large-scale struetures 
in the Universe. 

In 1992, the COBE (COsmic Background Explorer) satellite [5], [6] observed the 
CMBR with unprecedented precision and revealed for the first time that the level of 



S 2 = 




w 2 = 



E ab E ab + H ab H' 
6H 4 



(1) 
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the CMBR temperature fluctuations on large scales is as small as $f ~ 1(T 5 [7, 8]. 
After COBE many other ground and balloon born experiments [9], with higher angular 
resolution, confirmed this result and allowed us to probe the level of the anisotropies on 
a large ränge of scales. 

On large angular scales, the CMBR anisotropies (^-), are dominated by Sachs- 
Wolfe effect. This phenomenon, already deduced theoretically by Sachs & Wolfe [10], 
was used to compute the first-order perturbations in a FLRW universe with a flat 3- 
space filled either with dust or radiation. This is just one of the various possible sources 
of anisotropy, which occurs when there are inhomegeneities in the distribution of matter 
on the surface of the last scattering, that may produce anisotropies by the redshift or 
blueshift of photons. In this paper we compute the Sachs- Wolfe effect [10] for some 
anisotropic but homogeneous models (Kantowski-Sachs and Bianchi type-III) and find 
that under the assumption H a ~ these models give rise to the same Sachs- Wolfe effect 
obtained for FLRW universes. This is an interesting result witch teils us that CMBR 
observations on large angular scales will not be able to distinguish these anisotropic 
models from FLRW ones. 



2. The method 



As Collins and Hawking [11] pointed out, the number of cosmological Solutions which 
demonstrate exact isotropy well after the Big Bang origin of the Universe is a small 
fraction of the set of allowable Solutions to the Einstein equations. It is therefore prudent 
to take seriously the possibility that the Universe is expanding anisotropically and to 
investigate what effect anisotropic expansion will have on the angular distribution of 
background radiation [9]. In this work we show that, for large angular scales (i? > 10°), 
there exist homogeneous but anisotropic models, where the photons traveling to an 
observer from the last scattering surface encounter metric perturbations which cause 
them to change frequency, just like in the case of FLRW models. 

The metrics we consider are the Kantowski-Sachs and Bianchi type-III, given by 

ds 2 = -dt 2 + a 2 (t)dr 2 + b 2 (t)(d6 2 + f 2 (6)d<f) 2 ), (2) 

where 

, , [ sinö for Kantowski-Sachs 
| sinh 9 for Bianchi type-III 

We evaluate the Sachs- Wolfe effect [10, 12], assuming small perturbations in the previous 
metrics, and then integrating the geodesic equations for the CMBR photons along their 
paths, from the Last Scattering Surface (LSS) to the observer. In this work we account 
for the "kinematics effects" undergone by the free propagating radiation from the last 
scattering, in a perturbed universe, and for the "intrinsic effects" originated by the set 
of physical and microphysical processes related to the density perturbations in the LSS. 
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For simplicity, it is common to perform a conformal transformation § of the previous 
metrics and to work with the following metric forms 

d f = -dr ] 2 + dr 2 + ^l(d9 2 + f 2 (6)d(j) 2 ) (3) 
a 2 (r]) 

such that ds 2 : ds 2 = a 2 (r])ds 2 , since the null geodesics are preserved by this 
transformation. Afterwards, the results are transported to ds 2 metric. The metric 
ds is perturbed in the following way 

b 2 (n) r 

d f = - (1 + h m )dr, 2 + (1 + h u )dr 2 + -^v (1 + h 22 )d6 2 

+(1 + h^)f{6)d(t> 2 } - (V + h 10 )drjdr - ^p\(h 02 + h 20 )dr)dö 
J a(rj) 

0(77) 0(77) 

+ M/( ö )(/ ll3 + h^drd^+^Af^)^ + ^^rf^, (4) 
0(77) a 2 (r]) 

where /i a & are functions of time and position and such that /i a ;, <C 1. 
Considering the geodesic equation for the photons 

where ?7 a represents the photon 4-vector velocity components and u> the affine parameter 
associated to its trajectory, we may calculate this 4-velocity integrating the previous 
equation 

U a = ^Jg bc , a U b U c dw + (0) U a . (6) 

The term (o)U a represents the non perturbed photon 4-velocity components in the 
covariant form. 

A material observer, moving with some 3-velocity V, in a perturbed universe, has 
a 4-velocity given by (see Appendix A) 

V a ^(l- l -h m ,V\V 2 ,V^. (7) 

This observer measures a photon energy E 1 proportional to U a V a , 

E 7 oc Q f g bc , a U b U c dw + (0) C/ a ) (l - h 00 , . (8) 

Decomposing the 4-velocity in Taylor series and conserving only the first order terms 
U a ~ (o)U a + (i)U a (the calculation of (o)U a is in the Appendix B) and neglecting terms 
like h 00 (i)U a (because they are second order terms), the expression becomes 



Ey0C \l 9bc,oU b U c dw - -h 00 J gbc,o(o)U\ )U c dw + {0) U 

- \ { o)U h 00 + U-V+ l -J g bcm U\ 0) U c dwV l . (9) 

§ Owing to this transformation, r] is usually called the conformal time, and it is related to cosmic time 
by dt 2 — drj 2 a 2 {rj). 
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The last term of right hand side vanishes, because for i = r and i = =>- g~bc,i = and 
also because = 0. The term 

^oo / 9bcfi{ü) ub {o)U c dw = -hm j g bcfi{0) U'\ 0) U c dw 

= \hoo J g ee , v ( ( o)U e ) 2 dw = i/ioo ^ - ^ ((o)U 9 ) 2 dw, 

may be numerically computed. By Convention we use (') = ■£-. Uwe choose accurately 
the values of density parameters Vl M + (see Aguiar & Crawford [2]) Q M + — 1 
this integral may be neglected, because it is a second order term (see Appendix C). The 
first term in the right hand side may be decomposed by the following way 

\ J g bcß U b U c dw 

= \J 9bc,o(o)U b ( )U c dw + 9bc,o(o)U\i)U c dw 9bcß(i)U\ )U c dw 
= \\ 9bc,o(o) u \o)U c dw + J g gdv{ Q)U e {1) U ö dw. 



But 



/ heMo)U\i)U d dw = 2j ^ (± - ^) ( o)f/ e (i)f/ e ^, 



may also be neglected because the reason pointed previously. Thus, 

1 r 1 -> - 

Ejoc-J g b c,o(o) ub (o) UCdw + (o) u o ~ ^(o)U h 00 + U-V. (10) 

Calculating the ratio of energies in the emission instant (e) and in the reception instant 
fr) we obtain 



K (Wie -l + lh 00 \ e +(Ü-V) 



E lr (U a V a )\ r _ x + i hQolr +\}g bcm U\ ü) U c dw' 



(the symbol X\ A means that X is being evaluated at point A). Considering the 
approximation (l + X)^ 1 ~ 1 — X we have 



1 + \ NE +[V'V] r e + \J 9bc,o ( o)U b {0) U c dw, (12) 



e 



where [X] A = X(B) — X(A). The obtained results in ds 2 may be transported to ds 2 
using the relation E(t) = -^E{w), as we see in Appendix D. The photons redshifted 
from the last scattering (in e) until being observed (in r) may be calculated by the ratio 
of measured energies in emission and reception, 

\ e E(t r ) a e E(w r ) v 1 

because the wavelength of the photons is inversely proportional to its energy. On 
the other side the redshift is equal to the ratio between the black body associated 
temperatures in the emission and reception instants, 

P = * + l- (14) 
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Relating the previous equation and the Equation (13) and considering newly a linear 
approximation we obtain 



T r . ~ — T P 



a r 



1 " 2 M 



u-v 



9bc,o(o)U\ )U c dw 



T e 1 + 



TT 



therefore, 



^ 1 r/, 
TT " "2 [ho ° 



U-V 



1 r 

2 / ^ bc '° 



{0) U\ 0) U c dw. 



(15) 



(16) 



Now it's time to dehne the perturbations h a b- These fluctuations are gauge dependent. 
This correspond to the hypersurfaces choosing where these fluctuations are defined 
[13, 14]. We will choose a Newtonian gauge to allow a intuitive understanding. 
Considering only scalar perturbations the non zero quantities are, 

h 00 = 2^; h u = h 22 = h 33 = 2$. (17) 

As we have stated, these scalar quantities are function of time and position, \1/ = ^(rj, x l ), 
$ = $(77, x l ) and may be interpreted as Newtonian potential and a spatial curvature 
perturbation potential, respectively [13, 14]. 

For our models, with a proper choice of density parameters, these universes are 
approximately flat (Q + Üa — 1) and one may show that if one neglect the pressure 
(p = 0) one gets ^ = — $. Writing (16) in the Newtonian gauge we obtain, 

5T r 



u-v 



r\2 



8\2 



dw 



-Vo((o)C^) 2 + /iii,o((o)tr) 
b 2 (b 



The last term of right hand side may be neglected, as we show in Appendix C. If we 
consider vanishing pressure (\1/ = — $) we get 
5T r 
~Tr~ 



- m: - 



u-v 



or yet, 



5T r 

"TT 



U-V 



+ 2 



drj 



dw. 



e\2 



>dw, 



(19) 



(20) 



Note that the terms in square brackets of Equation (19)are equal to 1 as showed 
in Appendix B. Now, we should spell out the physical Interpretation of each one of 
three factors of the right hand side of previous equation. When matter and radiation 
decoupled, free CMBR photons, climbing the gravitational potential generated by 
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density perturbations, undergo a gravitational redshift, with corresponding loss of 
energy. The photon energy Variation in this process is given by the term [ty] r = 
(r) —ty (e) . The second term, Ü • V = Ü • (v(r) — V(e)j , corresponds to the Doppler 
effect induced by the relative motion of the observer in the emission and reception events. 
The last term teils us that the perturbing potential may vary between the emission and 
reception instants. 

The Doppler term has an observational meaning of a dipolar anisotropy on CMBR 
temperature and is usually removed from the equation to be treated aside. The last 
two terms are usually called the Sachs-Wolfe effect or also the integrated Sachs-Wolfe 
effect. For flat models of FLRW without cosmological constant, \1/ is time constant 
[12, 13], so the last term in the right han side of Equation (20) vanishes. In our case 
the cosmological constant is not vanishing, and indeed A will play an important role in 
our analysis. We consider values such that Qq + Qa — 1, taking into account recent 
observations [15, 16] which suggest Q ~ 0.3 and fl\ ~ 0.7. 

Equation (20) does not contain all physical processes which may generate 
fluctuations in CMBR temperature. It only accounts for kinematical effects undergone 
by the photons during their free propagation in a perturbed universe. So, we must also 
account for the intrinsic temperature fluctuations originated by the set of physical 
and microphysical processes, associated to density perturbations in LSS. Despite its 
youth, the Universe is already highly isotropic (shear a ~ 0). Then, for simplicity, we 
assume in this section that the Universe might be characterized by a flat FLRW model. 
Because the density fluctuations are very small, we may treat them in the context of 
linear theory of perturbations using Stefan-Boltzmann law 

p 7 = a SB T\ (21) 

where o~sb is a constant and p 7 is the radiation density. Differentiating this equation 
one easily obtain 

ST e 1 5o^ 

T e 4 p 7 1 ; 

At this time, when the Universe is very young, the total energy density is not only due 
to radiation. The baryonic matter play an identical role, and so the matter density p m 
is related with p 7 by 

^ - l S -Bl = 0, (23) 

Pm 4 p 7 

if the perturbation mode is adiabatic and on scales larger than the horizon at this 
time [17]. Nevertheless, if the perturbations are isothermal, the matter distribution is 
perturbed without making significant changes in the radiation density [17]. In this case 
and while the perturbations remain outward the horizon, we have 

^ ~ 0, (24) 
p 7 
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so, the temperature inside a perturbed region with a dimension greater than the horizon 
remain with constant temperature (5T e ~ 0). In summary, for perturbations on scales 
greater than the horizon we may write 

5T e { \ 5 -^- <= adiabatic pert. , . 

— = \ * Pm * _ (25) 

T e [0 <= isothermal pert. 

From the last expression we see that for adiabatic perturbations, the over-density 
(under-density) regions are intrinsically hotter (colder) than the LSS mean temperature. 
According with [12, 13], Sp m /p m = —2^ + Ö[(k/H) 2 ], where k is the momentum 
associated to perturbation scale and H the Hubble parameter||. The larger is the scale, 
the smaller is k, so, for perturbations greater than the horizon k <C H, the over-density 
locals coincide with the potential well, because, 
5T e 2 

jr * -gtte- (26) 

The measured temperature in a LSS point T e = T e (rj, x l ) may be related with the mean 
temperature < T e >, 

T e =< T e >[\ + 5 ^. (27) 

Using the expression T r = ^T e (l + ^p 1 ) an d substituting here the previous equation 
and retaining the first terms only, we obtain 

T, = ^<T e >(l + ^ + ^), (28 ) 

then, the total observed temperature fluctuation in r is 
5T 5T r 5T e 

~T ~ ~% + ~f7' ^ ' 

Substituting Equations (20) and (26) in (29) we have 

— = -y e -U-(V r -V e ) + 2 j —dw, (30) 

e 

where without loss of generality we put ty r — 0. The previous equation is valid if the 
Observation is made for regions with angular scales containing the horizon (i? > 3°) in 
recombination epoch. If the perturbations are isothermal, the temperature fluctuations 
coincide with ^ and are given by (20). 

Note that it is not necessary to include the second term of right hand side of 
Equation (30), since the dipolar anisotropy of the observer's motion is usually removed 
from the Observation data. Finally we obtain 

e 

which is the same expression obtained for FLRW universes for the same order of 
approximation, and for adiabatic initial conditions. 

|| Nevertheless our models have two Hubble parameters H a and Hb, they have practically the same 
value H a ~ H b ~ H due to our parameters choose. 
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3. Concluding Remarques 

We stress once more we should bear in mind that the assumption H a ~ Hb does not 
imply, by itself, an isotropic or even an almost isotropic metric, as is expressed by the 
growth of Weyl term in Equation (1), when we go back in time to the last scattering 
epoch. Although the £ 2 term remains at a low value from the present (Ejj ~ 0) to the 
last scattering epoch (£ 2 S ~3x 1CT 13 ), the Weyl term, as we computed in Appendix E, 
grows from Wq — 1CT 19 at the present to W? s ~ 1CT 1 over the same period of time. This 
shows the anisotropic character of these models in the past. Even though we impose a 
high level of isotropy at present time, its anisotropic behavior comes forward as we go 
back in time. Nevertheless, the growth of W 2 term does not affect decisively the first 
order computation of 5T/T term. 

Because the obtained expression (for Sachs- Wolfe effect) is the same as the one 
given for FLRW models, we may conclude that, these anisotropic models are also good 
candidates to the description of observed Universe provided we may assume H a ~ Hb 
and a particular choice of the density parameters: Q + J1a — 1) from the last scattering 
to the present, see Appendix C). This is another step taken in the same direction as in 
[2]. This is also in agreement with a known result: it is not possible to distinguish a 
Kantowski-Sachs model from the FLRW models, with the classical tests of Cosmology, if 
the Hubble parameters along the orthogonal directions are assumed to be approximately 
equal [1]. 

There are now many CMBR experiments in preparation that will allow to make 
much higher detailed observations than presently. Experiments of particular importance 
are satellites (MAP [18], PLANK [19]), and interferometers (AMIBA [20], AMI [21], 
CBI [22]). The first, have the great advantage of mapping the sky globally, but 
interferometers can achieve higher resolution and therefore probe to angular power 
spectrum to very high t. PLANK satellite best resolution is 5 arc minutes and present 
interferometers can now reach the 1 arc minute scale. 

In conclusion, Observation of Sachs-Wolfe effect plateau does not permit to 
distinguish between FLRW models and these anisotropic ones. To investigate this in 
more detail, it is necessary to consider and process the data from MAXIMA [23] and 
BOOMERANG [24] projects to regions smaller than the horizon at the last scattering 
(£ > 100, < 1°). Within this region of multipoles, perturbations are model dependent. 
Only with this Information we may conclude finally whether our Universe goes through 
or not by a real anisotropic phase. This will be the purpose of further work. 
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Appendix A. Estimate of an observer 4-velocity in a perturbed universe 

We will consider a material observer which is moving along Iiis worldline with a 4-velocity 
components 

dx a 



V a = 



dr ' 



(A.l) 



where r represents the proper time. Because he(she) is like a material particle, then he 
(she) is under condition V a V a = — 1. So, considering the absolute value we obtain 



9a b V a V b 



or also 



-(1 + h 00 )(V ) 2 + (1 + M^ 1 ) 2 + (1 + h 22 )-AV 2 f 



(A.2) 



+ (i + h^y-p mvy - (h 01 + h 10 )V u V L - -(h 02 + h 20 )V°V 2 
er a 

- -f(0)(h 03 + h 30 ) V°V 3 + -(h 12 + h 21 )V l V 2 + -f(9)(h 13 + h 31 ) VW 3 

(X (X CL 

b 2 

■ ' ' ' 1. (A.3) 



+^f(0)(h 23 + h 32 )V 2 V 3 



Making a first order approximation and since V 1 <C 1 and h^ <C 1 (with i = 1, 2, 3 and 
a, b — 0, 1, 2, 3), we may neglect products involving h a b and V 1 and also (V 1 ) 2 , resulting 



then 



'l + hoo)(V 



V°~l- -h 00 



0\2 



1, 



(A.4) 



(A.5) 



thus, the 4-velocity vector of an observer in a perturbed universe may be written in first 
order approximation by 



v a ~ {i-hi^v\v 2 ^ 



(A.6) 



Appendix B. Determination of photons 4-velocity in a non-perturbed 
universe 



The photons 4-velocity in a non-perturbed universe (o)U a may be reached recurring to 
geodesics equation in the metric ds 

<h)Üa _ 1 ~ fjb t-jc 

For (j) co ordinate we get 



(B.l) 



d (P) U <i> _ 1 ~ frb frc _ n 



(B.2) 
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because g bc ^ = 0, so (o)^ = const. =>■ ^ sin 2 9 (o)Ü^ = const. For 9 = =>- const. = 
=>- (o)^ — 0- For 9 coordinate we get % 

^ = lhc,e(o)Ü\o)Ü c = ^sinöcosö( (0) ^) 2 = 0, (B.3) 

because ( )C/^ = 0. Then, (0 )C/e = const. = v 7 /? =>• (o)^ e = a 2 /b 2 y/]3. 
For r coordinate we get 

^ = ^r(0)^(0)^ = 0, (B.4) 

because g bcr = 0, so (o)Ü r = const. = yfä =>- (o)Ü r = y/a. For a photon 
(o)Ü a{0 )Ü a = ( ( o)f/") 2 = ( (0) t/ r ) 2 + & 2 /a 2 ( ( o)^) 2 =» (o)^ = Jc7+aJ/¥ß. Thus 
the components Ü a of 4-vector are the following (^J a + a 2 /b 2 ß, y/a, a 2 /b 2 y/ß, o) , and 
C/ a will have the components (-^a + a 2 /& 2 /3, v^, v^, 0) . 

Determination of photons trajectory in a non-perturbed 
universe 

Let us consider that at the emission instant the photon has coordinates 
[x\ = t)e, x E , x 2 E , x|, w — 0) and at the reception instant has 

(x° E = t)r, x R — 0, x 2 R — 0, x R — 0, w — rjR — t)e)- In a general way we may consider 
the integral 

w 

x° = J {0) Ü°dw' + A, (B.5) 
o 

where A is an Integration constant that will be fixed by initial conditions of the problem. 
So, we get from previous equation that for the emission instant (w — 0), A — tje- Adding 
the reception data and substituting the (o)U° value, the integral becomes 

Vr-Ve I " — 

Vr-Ve= j Ja + ^ßdw. (B.6) 
o 

Düring the photons travel from LSS the Hubble parameters preserve the relation 
ä/a ~ b/b, that is a oc b. Thus 



a 2 



\Ja + —ß = const. (B.7) 

But by Equation (B.6) the constant must be 1, and (o)^° = 1 and (q)Üq — — 1- 111 
conclusion, ^Ü a = (1, yfä, a 2 /b 2 y/ß, 0) and ^Ü a = (—1, y/ä, y/ß, 0). 

f For Bianchi type-III, we will have sinhöcoshö instead of sin cos 0. 
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Appendix C. Numerical computation of an integral 

The Einstein equations for Kantowski-Sachs and Bianchi type-III metrics with perfect 
fluid, vanishing pressure and cosmological constant A, are 

« b b 2 kc 2 . 

2 1 + (C2) 

° + j+°j = A, (C.3) 
a o ab 

where p represents the baryonic matter density, G the gravitational constant and k 
takes on the values ±1 for Kantowski-Sachs or Bianchi type-III, respectively. The first 
integral of Equation (C.2) leads to 

^_Mi_k_ A 

b 2 - b 3 b 2 + 3 > y LA > 

where M 1 is an integration constant. We define the Hubble parameters corresponding 
to scale functions a(i]) and 6(77) as 

H a = ä/a e H b = b/b, (C.5) 

which may be used to define the density parameters, in analogy with which it is usely 
done in the FLRW models, 



b^H 2 ' b 2 H 2 3H 2 

From (C.4) we get a conservation equation 

n M + n k + n A = i. (c.7) 

Substituting (C.4) in (C.l) we get another relation that can be expressed, also, like a 
conservation condition 

% - Ü M + 2ft A = 2^, (C.8) 

where f2 p = M p /(ab 2 H 2 ), being M p a constant proportional to matter of the Universe 
(for details see Aguiar & Crawford [2]). The matter density parameter f2 may be 
obtained from the previous by 

q= _^p (c } 

l + 2#* K J 

Hb 

From equations (C.l), (C.4), (C.7) and (C.8) we also obtain the differential equation 

(CIO) 



dx %Mi-f) + Ovo(-i + V) + i^ 



dy ü Mo (l - y) + tt Ao (y 3 - y) + y 

where x = a/ao, y = b/bo and the under script denotes that the respective quantities 
are measured in actual epoch (see Aguiar & Crawford 
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Table Cl. Density parameters, relative difference between H a and Hb, and integral 
computation (J), for Kantowski-Sachs (KS) and Bianchi type-III (BIII) models. 





^M 




ü + ü Ao 




I/§< 






AH 

H a 




KS 


1 


< 1.7 x 10~ 8 


1 + 5.6 x 10" 


-09 


+5.9 x 10" 


-6 


-1.6 x 10" 


-6 


KS 


< 2 x KT 15 


1 


1 + 6.7 x 10" 


-16 


+6.8 x 10" 


-7 


-1.4 x 10" 


-6 


KS 


< 0.3 + 7.0 x 10~ 9 


0.7 


1 + 2.3 x 10" 


-09 


+3.5 x 10- 


-5 


-1.7 x 10" 


-6 


KS 


0.3 


< 0.7 + 7.0 x 10~ 9 


1 + 2.3 x 10" 


-09 


+3.5 x 10" 


-5 


-1.7 x 10' 


-6 


BIII 


1 - io- 10 


> 9.9 x 10" 11 


1 - 3.3 x 10" 


-13 


-1.4 x 10" 


-5 


+1.3 x 10- 


-6 


BIII 


> 9.8 x 10~ 14 


1 - 10~ 13 


1 - 6.7 x 10" 


-16 


-6.9 x 10" 


-7 


+1.3 x 10- 


-6 


BIII 


> 0.3- 10~ n 


0.7 


1 - 3.3 x 10" 


-12 


-1.1 x 10" 


-5 


+1.8 x 10- 


-6 


BIII 


0.3 


> 0.7- 10" 11 


1 - 3.3 x 10" 


-12 


-1.1 x 10" 


-5 


+1.8 x 10- 


-6 



We are assuming that the two Hubble parameters along orthogonal directions are 
presently indistinguishable (H ao = H bo ). With this restriction and using the previous 
equation we expect to compute numerically the integral 

From Appendix B we take (o)U 9 + and as dn/dw = (o)U v = 1 we get 

ä/a\ b 



i = ßj 




1 - 




6/6, 

The numerical Integration is made between 1/1000 to 1. 

If we choose the density parameters Qm + ^a () very near the unity, (under or 
over one), we get low values for the integral I, for Kantowski-Sachs or Bianchi type- 
III, respectively. In the box below we can see the numerical Integration result for 
several combinations of density parameters Qm an d ^a ; without the multiplicative 
constant al/b^. We computed also the percentual difference between Hubble parameters 
at z = 1000 (AH/H a = (H a - H b )/H a ). 

With the obtained values, if üq ~ 6o, the integral / is second order, then we neglect 
it, so the quantities ^hool and 2^U 9 I must also be neglected. 

+ The ß value may vary between to 1 and teils us about the velocity percentage that a particle has in 
an angular direction (in this case the photon). Thus, if the velocity is purely radial ß = 0, by Opposition 
when the velocity is purely angular we have ß = 1. 
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Appendix D. Relations between 4-vectors expressed in a conformal and 
non-conformal metrics 

Let us consider the metrics ds 2 and ds 2 (equations (2) and (3)). As we have ds 2 = a 2 ds 2 , 
then 

#00 = 9oo and g ü = a 2 g iU 

with i = 1,2,3. On the other hand, the conformal parameter obeys to dt = adrj. 
This allows us to establish the relation between the photon 4-vectors for these frames. 
Labeling x a (£) and x a (w) the photon coordinates for these two frames, then the 
respective 4-vectors will be 

U a = ^- and Ü a = ^- (D.l) 



thus, 



and 



d£ dw 

U° = — = — (D.2) 

di { } 

= dx^ = dr]_ = l(tid£_ = 1 d£ , ß ^ 

dw dw a d£ dw a dw ' 



and the covariant component is 



U = 9o U° = g 00 -^U° = -^U . (D.4) 

aaw a dw 



As x l ([£) = x l (w) and U % = it comes to 



and 



Thus, 



dw dt, dw dw 

Üi = guÜ 1 = \gu^U l = \p-Ui. (D.6) 
er dw a z dw 

Ü a Ü a = Ü°Ü + Ü% = -^U°-^U + ^-ipL^-Ui 

a dw a dw dw er dw 

=H0 (u ° Uo+u ' u ' )= H& u ° u °- (D - 7) 

Although the photon 4-vectors must have null norms in the two frames, its ratio is 
defined by 

üw a -rf[fo) ■ (D - 8) 

It is easily seen that the temporal components*, in the conformal and expansion 
frames, we should have 

t/° = ^ = l and U° = ^ = -, 
dw ez£ a 

* Note that using the geodesic equation = \gbc.aU h U c , we obtain U r = y/ce, Ug — \fß and as 



U a U a = we get U a = (1/a, yß/a 2 , y/ß/b 2 , 0) and U a = (-l/a,y/ä,y/fi,0), because Ja + ^ß = 1. 
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thus, the i] and w parameters are both conformal, and because this frame has no 
expansion, the photon 4-vector temporal component (Ü°) has no changes, in Opposition, 
in the expansion frame this component decrease in an expansion phase. From the below 
relations and from (D.3) we obtain 

^ = a 2 . (D.9) 

dw 

Let us now obtain the relationship between the 4-velocity vectors in these two 
frames, for material particles. 

Let it be V a and V a the respective 4-vectors for the expansion and conformal 
frames. For material particles, these 4-vectors have the same scalar product, which is 
then frame independent 

V a V a = -1 = V a V a , (D.10) 

thus, 

goo(V ) 2 +9u{V i f = g 00 (V ) 2 +g ü (V 1 ) 2 - (D.11) 

As goo = #oo, 

9oo(Vy + g ü (Vy = g 00 (V ) 2 + - 2 g ü {V'f. (D.12) 

er 

The relation below allows us to conclude that if 

V° = V° then V 1 = -V\ (D.13) 

a 

Indeed, as we know 

dx° = adx° and dx i = dx 1 . (D.14) 
Then, we may write 

V a = a ^- and V a = ^-, (D.15) 

where T\ and r 2 are, respectively, the proper time for the material particles in the 
considered frames. Using relation (D.10) one gets 

-dr\ = -(dx ) 2 + a 2 {dx 1 ) 2 and - dr 2 = -(dx ) 2 + a 2 {dx 1 ) 2 . (D.16) 

By Equation (D.14) and calculating the ratio of T\ and r 2 we obtain 

dr[ -a 2 (dx ) 2 + a 2 (dx 1 ) 2 
d^l ~ -(dx ) 2 + (dx 1 ) 2 

or also 



(D.17) 



dn = adr 2 , (D.18) 

the choosing of the positive sign is obvious. Considering the first equation of (D.14), we 
may write 

dx° dn _ jlx®_ 
dn dr 2 dr 2 ' 
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or 

V°a = aV° =>■ V° = V° (D.20) 
as we claimed above. So, 

V = V and K = gaV 1 = gacf-V* = aVt. (D.21) 

a 

An observer measures the photon energy for these two frames as the time component 
of the photon 4-vector in his proper frame, 

E(t) oc U a V a = U°V + U% 

= -Ü°V + \ü l aVi 
a er 

= -Ü a V a oc -E(w), (D.22) 
a a 

as the proportional constant is the same, then 

E(t) = -E(w). (D.23) 
a 

Appendix E. Computation of £ 2 and W 2 

For Kantowski-Sachs and Bianchi type-III models, the shear tensor is such that 

a ab a ab = 2a 2 = 2 -{H b -H a ) 2 (E.l) 

as usually, H a and H b are the Hubble parameters in orthogonal directions. So, taking 
the definition given in Equation (1) 

S 2 1 ( H b -H a 
9 V H a 

because H ~ # a ~ if b . As we see in Table Cl (Appendix C) (^) ~ 1.7 x 10~ 6 for 

last scattering and (^7) ~ presently. Thus, S 2 S ~ 3 x 1(T 13 and S 2 ~ 0. 

Given a 4-velocity field, £7 a , of an observer set, one may split the Weyl tensor C abcd 
in the electric (E ab ) and magnetic (H ab ) parts, 

E ac = C abcd U b U d 1 (E.3) 
H ac = \c: bcd U b U d , (E.4) 

where C* bcd is the dual of C abc( i. For a comoving observer, the magnetic part is null and 
the electric part is such that 

2 

n ^ nh 1 / a ab b k b 2 \ , 



(E.2) 



Using Equations (C.2) and (C.3) we get 



= !(;-!)'• (K6 > 
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redefining the scale factors as x = a/a and y = b/b and using Equations (2.12) 
and (2.22) of Aguiar & Crawford [2] we derived x and y and after some algebraic 
manipulation we get 

(E.7) 



rp rpab Tri 

h ab h - -H ho 



n \ 3 1 i i oQ 1 2 Ha ° 1 

"M —7 ö + ^"A ö - ^77 ö 



E ab E ab = -<(n M „ + ^A - l) 2 - (E.8) 



For to, x — y — 1 and as if ao = _£4 we have 

2 

3^ 

We chose Qjw and Qa Parameters such that its sum was near the unit, in order to have 
high a level of isotropy 

I^Mo+^Ao-ll^lO- 9 . (E.9) 

Thus, for to we have 

w 2 = ^(n Mo + n Ao - 1) 2 ~ io- 19 , (e.io) 

and for the last scattering epoch x ~ y ~ 1/1000 we get 
1 / + 
9 V ÜF9 
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